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Abstract 

While the linearity of the Schrodinger equation and the superposition prin- 
ciple are fundamental to quantum mechanics, so are the backaction of mea- 
surements and the resulting nonlinearity. It is remarkable, therefore, that 
the wave-equation of systems in continuous interaction with some reservoir, 
which may be a measuring device, can be cast into a linear form, even after 
the degrees of freedom of the reservoir have been eliminated. The super- 
position principle still holds for the stochastic wave-function of the observed 
system, and exact analytical solutions are possible in sufficiently simple cases. 
We discuss here the coupling to Markovian reservoirs appropriate for homo- 
dyne, heterodyne, and photon counting measurements. For these we present 
a derivation of the linear stochastic wave-equation from first principles and 
analyze its physical content. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



Perhaps the most fundamental principle underlying quantum mechanics is the superposi- 
tion principle for quantum amplitudes which is equivalent to the linearity of the Schrodinger 
equation. However, as is well known, the time-evolution of the state vector of any sys- 
tem is governed by the linear Schrodinger equation only as long as the system is decoupled 
from its environment, and, in particular, from any extraneous measuring apparatus. E. g. if 
during a short time interval such an extraneous coupling is introduced, which performs a 
measurement of an observable X without destroying the system, the linear time-evolution 
is interrupted; the wave-function is replaced by a new wave- function chosen at random 
among the eigenfunctions \ip n ) satisfying X\ip n ) = X n \ip n ) with the eigenvalues X n ; the 
probability for \ip n ) is given by 

Pn=\(^nM 2 , 716 1 (1.1) 

where I is the set of all possible quantum numbers, here taken as discrete, for simplicity. 
This step in the time-evolution is therefore random and nonlinear. In order to formulate 
it explicitely let us introduce a classical random variable N, which only takes the values 
iV = n, n 6 I with probabilities 

P(N = n) = p n , $> n = l. (1-2) 

Then we may write for a measurement 

IV) -> W) = \M. (i.3) 

If similar measurements are performed repeatedly in time, say at the discrete time ti, then 
the linear deterministic time-evolution between successive measurements is repeatedly inter- 
rupted by steps like ( |1.3| ) with more random variables iVj whose probabilities follow again 
from eq. ( |1 . 1| ) and are therefore correlated, in general. 

The possibility to describe measured systems, or more generally systems coupled to their 
environment, by stochastic wave-functions similar to (|1.3|) is of course well known. By 
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the work of many authors over the years this possibility has been turned into an efficient 
tool, useful in particular for numerical simulations, by deriving equations of motion for 
the wave-functions of systems continuously coupled to their environment. All of the early 
and some of the more recent work in this area is purely phenomenological (see e.g. 
and often motivated by the desire to find generalizations of the Schrodinger equation on a 
fundamental level including wave-function collapse. Some of the more recent work, however, 
stresses microscopic derivations, in certain weak coupling and Markovian limits, from the 
Schrodinger equation for the coupled system and reservoir (see e.g. [|]-f|). As one may 
expect from eq. (|1.3| ), the resulting wave-equations generalize the Schrodinger equation for 
isolated systems in two ways: they contain in an explicit way classical stochastic variables 



analogous to N in eq. (|L3"1), and they are nonlinear, analogous to the nonlinear dependence 
of N on 1*0) in eq. (|L3]) via its distribution Q, (|i~T| ). 

The appearance of stochastic elements in the generalized Schrodinger equation may look 
unusual at first sight but should not cause surprise, as it corresponds directly to the observed 
stochastic behaviour e. g. of measured quantum systems. However, the appearance of a 
nonlinearity in quantum measurements, even though well known to all physicists, may still 
be considered as alarming, because it seems to change the foundations of quantum mechanics: 
In the real world, perhaps with the exception of the universe as a whole, there aren't any 
systems which are completely decoupled from their environment. This seems to lead to 
the alarming conclusion that linearity and the superposition principle are only approximate 
concepts for any system in the real world. Indeed, any finite real-world system in an excited 
state will sooner or later decay via one or several random nonlinear steps similar to eq. ( p..3|) . 
But if the nonlinearity of the environment-coupled reduced time-evolution must be accepted, 
the role of the superposition principle becomes obscure. E. g. it may no longer be possible 
to analyze the fate of an initially given linear superposition of states by studying the fate of 
each of its constituents separately. 

In view of this fundamental problem raised by the nonlinear stochastic evolution ( |1.3| ) 
it is very remarkable that the stochastic nonlinear wave-equations, describing a broad class 



of systems in continuous interaction with their environment, can actually be cast into a 
linear form, containing the same information. These linear equations still contain classical 
stochastic variables, however now distributed independently of the wave-function. The noise 
they describe is therefore independent of the system under study and may be considered a 
property of the environment alone. 

The existence of these linear versions of the stochastic wave-equations has repeatedly 
been noted in the literature [jsHT^I, however their fundamental relevance in connection with 
the superposition principle has, to our knowledge, not been stressed. In the present paper 
we wish to derive and study in a unified manner the linear stochastic wave-equations corre- 
sponding to homodyne, heterodyne and photon-count measurements in quantum optics. 

These measurements have been considered before by Wiseman and Milburn |||7| using 
nonlinear stochastic wave-equations. A microscopic derivation of a linear stochastic wave- 
equation corresponding to homodyne measurements was presented by Belavkin [|]|§. Our 
derivation starts from the same assumptions, but appears to be simpler, at least to us. 
Photon-counting has been discussed in the framework of phenomenological nonlinear and 
a phenomenological linear stochastic wave-equations in papers by Barchielli |T2|,|TSJ and by 
Belavkin |1U|| . The linear stochastic wave-equation for photon-counting, which we shall 
derive in the present paper is different from the phenomenological version used in their 
work. Stochastic Schrodinger equations with quantized noise terms have been discussed in 
detail in the work of Gardiner et al [|14j]. While these authors consider equations in the large 
Hilbert space of system and environment, our work here is concerned with wave-equations 
in the Hilbert space of the system without environment. An extensive review of the use of 
nonlinear stochastic wave-equations in quantum optics and a large list of references is given 

mag. 

The derivation and use of stochastic wave-equations, both nonlinear and linear, is so far 
restricted to Markovian reservoirs and we shall also adopt this framework here. It remains 
an interesting open problem for future work to examine also non-Markovian extensions. 
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II. LINEAR STOCHASTIC DESCRIPTION OF A GENERAL QUANTUM 

MEASUREMENT 



Here we shall present the general idea to be applied in the subsequent sections for the 
case of a general quantum measurement. It is our aim to show how the stochastic description 
to be applied in subsequent sections is firmly rooted in the basic formalism established by 
von Neumann (1931) Q, Wigner (1963) f[7| and others. 

The nonlinear stochastic change of a wave-function due to the idealized instantaneous 
measurement of a quantum observable was recalled in eq. ( |1.3|) . The nonlinearity resides 
in the dependence of the probability p(N = n) = \{ip\ip n )\ 2 on \ip). We now turn to a 
stochastically equivalent linear description, i.e. one which gives the same results for physical 
probabilities. 

We may begin by noting that the transformation (|1.3|) would trivially be linear if the 
distribution of the random variable iV would be independent of the state \ip). Let us therefore 
introduce a new classical random variable N also distributed over all quantum numbers n e I 
with some fixed but arbitrary distribution p n satisfying 

p(N = n) = p n > 0, p„^0Vn, = 1. (2.1) 

nei 



Let us then consider the replacement of the nonlinear transformation (|1.3| ) by the linear one 

\1>)^\#n) = -^\M(M'I>) (2-2) 

where N is distributed according to Q2.1|) . We note that |^) is not normalized, 

fefe) = ^-K^}|Vl (2-3) 
Pn 

i.e. the transformation ( |2.2|) is not unitary. However, we have instead 

EPnWnK) = I- (2-4) 

The probability Px{X n ) to observe the eigenvalue X n of the measured observable X is fully 
encapsulated in the linear stochastic transform of via the formula 
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Px{X n )=p n ^' n W n ). (2.5) 

Using eq. ( |2.2j ) in eq. ( |2.5| ) it can be seen that the usual formula Px(X n ) = |(V ; ra|V ; )| 2 = p n 
results and, furthermore, that p n drops out completely and may remain arbitrary, as long as 
the conditions ( |2.1| ) are satisfied. E.g. one may replace p n by a constant and interpret the 
corresponding equidistribution of N over the set I as reflecting our complete ignorance of the 
outcome of the measurement prior to the specification of and prior to the measurement. 
Correspondingly, eq. Q2.5|) represents the transformation from the 'input' distribution p n , 
reflecting our ignorance, to the 'output '-distribution p n produced by the measurement. This 
particular way of thinking about eq. (|2.5| ) will be useful lateron, where descriptions analogous 
to eq. Q2.2|) are applied to various concrete measurements which are performed continuously 
in time. In fact, concrete expressions for the probability p^ are obtained in these concrete 
models. 



From eq. (|2.5| ) all moments of the measured observable may be calculated. Furthermore, 
it is straightforward to check that for measurements performed successively in time multi- 
time joint probabilities of the observables and the correlation functions associated to them 
may also be calculated by the repeated use of eq. ( |2.2| ) for each measurement. Hence, 
eq. ( p.2|) provides a description of the measurement of X which is just as complete as the 



nonlinear description ( |1.3|) , but linear. E.g. if the transformation (|2.2|) is known for two 



wave-functions \i/ja) and \iPb) then we have for their linear superposition 

C a \^a) + C B \ij B ) -> CU|Va> + CbWb) 

4=|^> (Ca{iPnM + Cb&M) ■ (2.6) 



'Pn 

We note that this expression still describes the stochastic 'collapse' of the superposition 
to a single state \ip N )- Yet the use of eq. ( |2.6| ) in the formula ( |2.5| ) correctly predicts the 
interferences of the superposition observed in the measurement of X based on the behavior 
of the individual components of the mixture. 

We may ask now whether the wave-function \if)'jj) is physical, i.e. whether it corresponds 
in any sense to the actual state of the reduced system. This question is difficult or perhaps 
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impossible to answer, even for the wave-function of closed systems, which may well be 
considered as being devoid of physical reality and no more than a mere computational tool 
to obtain physical probabilities and their moments via some definite algorithm. The same 
may be said for the wave-function IV^-), the only difference being that the usual algorithm 
for obtaining probabilities is slightly changed according to eq. fl2.5p . 

In order to derive eq. (|2.2|) from the Schrodinger equation we may use the method of von 
Neumann (1931) |16| and Wigner (1963) fl?]]- They show that an uncorrelated initial state 
|^) of the measured system and |0) of the measuring apparatus evolves according to 

I^>®I0>^E(^>I<M®I<M ( 2 - 7 ) 



into an entangled state. The entanglement in eq. ( |2.7|) is real and could in principle be 
observed if further joint measurements on system and apparatus are made. However, as 
long as further measurements are made only on the system or only on the apparatus the 
relative phases between the components of eq. ( |2.7| ) for different n are unobservable, because 
both the \ip n ) and the \4> n ) are orthogonal for different n, the \ip n ) because they correspond 
to different eigenvalues X n , the |0 n ) because they are macroscopically distinguishable states. 
Hence, as far and only as far as such restricted further measurements are concerned, a super- 
selection rule holds for n, (i.e. states of different n are not coherently superposed) and the 
right hand side of eq. (|2.7f ) is equivalent to a mixture 

£|(Vg^>| 2 |^n><V>n|® \<t>n)(<t>n\- (2.8) 
n 

This restriction of the type of further measurements allows to circumvent von Neumann's 
and Wigner's argument that a pure state cannot evolve into a mixture, an old stumbling 
block in quantum measurement theory. Tracing out the apparatus we end up with a mixed 
state, which is equivalent to the description ( |2.2| ) and ( |2.1| ) for arbitrary p n . 
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III. HOMODYNE MEASUREMENTS 



Let us briefly recall the principle of homodyne measurements. A quantum oscillator or 
mode, say the electromagnetic mode of a cavity at a given frequency but with (in comparison 
slowly varying) time-dependent amplitude, radiates into the vacuum, where the classical field 
of a local oscillator at the same frequency with stable phase and time-independent amplitude 
is superimposed and the beat intensity is detected. In balanced homodyning a difference 
scheme is employed to get rid of the dominating field intensity of the local oscillator in the 
output. 



The input-output formalism |18| is perfectly adapted to describe this situation: The 
interaction of the system, i.e. the quantum oscillator, with an input field generates an output 
field depending linearly on the amplitude of the oscillator via a coupling constant. In the 
case at hand the input field are just vacuum fluctuations entering the quantum oscillator by 
the time-reverse of the channel through which it radiates into the vacuum. For a differential 
time-step we may write 

dB out {t) = y/ja(t) + dB m (t). (3.1) 

Here a(t), o/(t) are the mode operators of the oscillator in the interaction picture (i.e. with 
the oscillator frequency lu splitt off); y^j is a coupling constant which we take positive by 
appropriately adjusting the origin of the phase of a; dB in (t) represents the future directed 
increments of the quantized input-field which is taken as quantized Gaussian white noise 
1 19|1 . Therefore 

[dB in {t),dB in {t')\ = = [dBi{t),dBl{t')\ 
[dB in (t),dB} n (t')\ = for t^t' 
[dB in (t),dB} n (t)] = dt 

[dB- m (t),a(t')} = [dB} n (t), a(t')\ = for t > t' . (3.2) 

Acting on the vacuum we have dB} n (t)dBi n (t) = 0. As the increments dB m (t), dB- m (t) 
are causally disconnected with the quantum oscillator at time t all their commutors with 



the oscillator variables at the same time vanish. The operator dB ont (t) in eq. ( |3.1| ) is the 
output field generated in dt by the interaction of the system oscillator with the vacuum. 
The Hamiltonian for this interaction is 

H int (t) dt = H (a, a\ t)dt + i^(adB\ n (i) - dB in (t)a r ). (3.3) 

Here H (a, a' , t) is the part of the Hamiltonian responsible for the generation of the quantum 
oscillation, in the interaction picture. With the unitary operator 



U int {t) 



exp (-ij Q dsH int (s] 



(3.4) 



where [...]+ denotes time-ordering, we may then write 

dB out (t) = uUt)dB in (t)U int (t) (3.5) 



which gives back eq. fl3.1|) . The interaction Hamiltonian and eq. (|3.1| ) are therefore indeed 
consistent. As dB in {t) and dB out (t) are connected by a unitary transformation, the commu- 
tation relations of dB out (t), dBl ut (t) at equal time coincide with those of dB- m (t), dB- m (t). 
The signal measured in the homodyne experiment may now be expressed in terms of the 
output field and its adjoint as 

\(3\&(t) = fdQ(s) 
Jo 

\[3\dO(t) = |/3| (dB out (t)e-^ + dBl t (t)e^) . (3.6) 

Here (3e lv is the complex amplitude of the classical local oscillator. If this is expressed in 
terms of the input field we find 

dG(t) = 2y/jX(t)dt + dE(t) (3.7) 

with the input field 

dE(t) = {dB in {t)e~ ilp + dB} n (t)e i(fi ) (3.8) 

and the system observable 



X(t) = -(ae-* p + a t e i *) (3.9) 

measured in this experiment. 

The commutation relations ( |3.2| ) imply that 

[dS(t),dS(f)] = (3.10) 

and it follows from eq. Q3.7D and eq. ( |3.2D that 

[dO(t),dE(t)\ = = [dO(t),X(t)\ . (3.11) 

Therefore, there is a representation in which all operators in eq. ( |3.7|) are simultaneously 
diagonal and may be replaced by their simultaneous eigenvalues, which we shall denote by 
corresponding lower case letters. Thus 

d0(t) = 2y/^x(t)dt + d£(t). (3.12) 

The measured output field is therefore generated by the incoming vacuum fluctuations and 
the instantaneous eigenvalue of the X-quadrature component of the quantum oscillator. The 
relation (|3.12|) implies that the output process d9(t), like the input process obeys the rule 
of Ito calculus (d9(t)) 2 = dt. 

Eq. (|3.12|) defines the output-process d6{t) conditioned on the input-process d£(t) and 
the value of x at time t. In particular, if the distributions of d£(t), p^(d^(t)), and of x at t 
are independent, as will turn out to be the case lateron, the distribution of d9(t) conditioned 
on x(t) is given by 

p e {d9{t)) = Pz{d6{t) - 2^x(t)dt) (3.13) 

It remains to derive the stochastic Schrodinger equation. The generator of the total 
time-evolution follows from 

tZ|V(t)> = dU int (tM(0)) = dU iQt (t)UUt)\m) (3-14) 

as 
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L(t)dt = dU int (t)Ul t (t) 



(3.15) 



where dUi n t(t) is calculated from Ui nt (t) by using the Ito calculus for the stochastic differential 
dB in (t) and its adjoint. The result is 



L(t)dt = —iHo(a,a\t)dt + ^(adB\ n {t) — dB in (t)a^) — ^-a)adt 
— — (a'a + aa')dBj n [t)dB- m (t). 



(3.16) 



Initially the total wave-function 1^(0)) ® |{0}) is the vacuum state |{0}) of the field. The 
total wave function at time t may be represented as the time-ordered product over discretized 
time 



H (1 + L(ti)dti 



U=0 



|V(0)>®|{0}>. 



(3.17) 



Then we may use the fact [dB- m (t), dB-^it')] = for t ^ t' to bring all dB(ti) to the right 
where they are annihilated when hitting |{0}). By the same token we may replace dB\ n (ti) 
by 



(3.18) 



which may be further replaced by d^{ti)e ltp in the ^-representation in which all oE(i) are 
simultaneously diagonal. Thus Lit) in eq. ( |3.17j ) is reduced to 



, 7 + _ • 

L(t)dt = —iHo(a, a' , t)dt a'adt + y/^ae tip d£,(t). 



(3.19) 



At this stage there are no operators left in eq. (|3.17|) which act on |{0}), which may therefore 
be divided out. The essential point here is that the component of the wave-function repre- 
senting the environment (i.e. the vacuum) is not eliminated by taking a trace, but rather 
by division. Therefore no information is lost in this step. 



For a single infinitesimal time-step the expression ( 3.17 ) now reduces to a linear stochastic 
wave-equation 



-iH (a, a\ t)dt — —a)adt + ^ae~ tip d^{t) 



(3.20) 
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In order to recover the full wave-equation, including the vacuum which was eliminated , both 
sides of eq. ( |3.20| ) have to be multiplied by <fo({£}), the vacuum |{0}) in the ^-representation. 
This functional is Gaussian and its absolute square 



(3.21) 



gives the classical Wiener measure of the white noise process £(s) with = fod£(s), 
(d£(s)) 2 = ds. We shall denote by }) the Wiener measure of £(s) for < s < t. It 

should be noted that \tp^(t)) is n °t normalized. However, the statistical operator defined 
by p(t) = J djjif {{^}) \ip^{t))(ip^{t)\ is correctly normalized. If one wishes the usual master 
equation for p(t) can be derived from eq. ( |3.20| ), but this is of no concern here. 

From the unitarity of the total time-evolution with H mt it follows that all expectation 
values of oscillator variables Y are given by 



<n*)> 



\w({0)\ 2 d{OiMt)\Y\Mt)) 



(3.22) 



More generally all time-ordered multi-time correlation functions of not necessarily hermitian 
oscillator variables Yi, Y 2 , Y$ etc. are given by 

(Ylfo)Y}(t 2 ) . . . Y r l(t' n )Y m (tm) ■ ■ ■ Y 2 (t 2 )Y 1 (t 1 )) 

= fdfi w m) (Mt)\YlKl(& ~ %)Y} ■ ■ ■ K\{t' n - Ci)^^(4 - t m )x (3.23) 

xY m K^t m - t m _i) . . . Y2K&2 - ti)Fi|^(ti)) 
where t' n > £^_ x > ... > t' x and t m > £ m _i > ... > t\. Here K^(t) is the (non-unitary) 
time-evolution generated by eq. (|3.20|) which may be written as 



K&) 



exp 



it -7 t 
H — t—a'a — 1- 



1a 2 e- 2iLp ) ds + yFjae-* I d£{s) 
2 / Jo 



(3.24) 



K^(t) is the adjoint of K^(t). The kernel Lit) in the middle of eq. ( |3.23|) is 



Ldt) 



Kz(t) if t > 
K\(t) if t < 
1 if t = 0. 



(3.25) 
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Thus, the full information about the reduced dynamics of the quantum oscillator is contained 
in the linear stochastic wave-equation ( |3.20| ). Even though K^(t) is nonunitary, it is unitary 



under the integral over g?/4^({£}), 

Jdfj?at})KlW&) = l- (3-26) 

This may be checked explicitely but also follows immediately from the unitarity of the 
time-evolution of the total system. 

The full information about the output process d9(t), regardless of the internal dynamics 
of the system, is also contained in eq. (|3.20|) . Indeed, the measure dfi t ({0}) defined by 



putting Yi = l = Y} in eq. (ET25D 



dfH({9}) = d^({0})(Mt)\Mt)) (3-27) 

can be identified with the measure of the output process for < s < t. We note that 
eq. ( |3.27| ) is just a variant of eq. ( |2.5| ) of section 2. To see the relation with the output 



process one computes the conditional measure 



in terms of the conditional Wiener measure 

p w mt)) = T.r/r^v (3-29) 

of the input process. One obtains 



dfiVm) 



(jnuw W tM n\\(-i , d (M t )\M t )) \ 



(3.30) 



From eq. ( ^.20| ) we obtain 



d(Mt)\Mt)} ^ (MQl^ + aWlMt)) ^) (3 3i) 



(Mt)\Mt)) (Mt)\Mt)) 

Recalling that p w (d$(t)) is just a Gaussian, formally written as 



p <*<*» = ^^(-drj- (3 ' 32) 
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we obtain equally formally, using the Ito calculus with (d6(t)) 2 = dt 



P(d9(t)) = ^== ex P 
V2ndt 



(3.33) 



Therefore 

d${t) = 2y/j(X(t)) e dt + d£(t) (3.34) 
where d£(t) is the Wiener increment and where we defined 

(xm (Mt)\x\Mt)) , . 

Eqs. (|3.34|) and ( |3.33| ) may be compared with eqs. fl3.12|) and Q3.13D , respectively. In the latter 
equations the output process d9(t) was conditioned on the simultaneous eigenvalue x(t), 
while in eqs. ( |3.33| ), (|3.34p it is conditioned, via (X(t))g, on the values of d6(s) for preceding 



times s < t. Therefore, because their probability distributions differ, the increments defined 
by eq. (|3.12| ) and (|3.34|) are not the same. (X(t))o is the quantum expectation value of the 



X-quadrature component for a given realization of the output process at times prior to t. 

Let us evaluate the multi-time correlation functions (d8(t n ) . . . d9(t\)) of the output pro- 
cess implied by eq. (|3.27|) and (|3.34j ). In the measure (|3.27|) of 6{t) we use the representation 



\Mtn)) = K e (tM{0)). (3.36) 

In (|3.27|) the integration is performed over the Wiener measure. In order to achieve this 
explicitely, for each time t{ appearing in the multi-time correlation functions the infinitesimal 
step dti in the functional ( |3.36| ) of dd(ti) is represented explicitely as (1 + L(ti))dU like in 



eq. ( ETT7D with 

L(ti) from eq. ( |3.19| ), making the dependence on d8(ti) explicit. The Wiener integrals 
over all the dB(tj) can then easily be performed with the result 

(dO(t n ) . . . d9(tx)) = (2^) n dt n . . . dt x (: X(t n ) . . . X(tx) :) (3.37) 

where : : denotes normal and time-ordering. The result ( |3.37| ) holds for arbitrary correlation 
functions as long as £j 7^ tj. For all pairs tj, tj with ti = tj one has to take d9(ti)d6(tj) = dti 
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summing over all possible pairings and consider the correspondingly reduced correlation 
function. 

In summary, the linear stochastic wave-equation has the following uses: 

(i) Solving it with d£(t) representing Wiener noise it can be used to obtain all correlation 
functions of system (i.e. oscillator) variables via eqs. (pC 



(ii) Solving it after replacing d£(t) — > d6(t) representing the yet unknown output noise 



(i.e. taking d6(t) arbitrary) one obtains the conditional quantum expectation (|3.35|) of 



the X-quadrature as a functional of the output noise from eq. ( |3.35| ). All correlation 



functions of the measured X-quadrature are then given by the correlation functions 

( CT of (H§. 

(hi) Finally the solution of (i) gives the complete measure of the output field when inserted 
in eq. fl3~27l) . 



If one is interested in generating stochastic numerical samples of (X(t))g and hence dO(t) 
one first makes the replacement d£(t) — > d6(t) in the linear stochastic wave-equation and 
then reexpresses the unknown output process d9(t) by the known Wiener process d£(t) via 
eq. Q3.34Q . However, this step generates a nonlinearity. The superposition principle is then 



lost in these numerical simulations. For nonlinear stochastic wave equations and numerical 
examples see [|,|H].[2(].(21] . 



IV. EXACTLY SOLVABLE EXAMPLES 

We consider two examples where eq. (|3.20| ) can be exactly solved. The first example 
is a harmonic oscillator periodically driven by some external field with amplitude F, and 
frequency u = ujq + 5 in the rotating wave approximation. Thus 

H (a,a\t) = 1 (F*e i5t a - Fe^V) . (4.1) 

With a coherent initial state 
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|^(0)) = |a), 



(4.2) 



where 



a\a) = ct\a), 



the solution of eq. ( p.2U| ), after some calculation, is obtained as 

1 



|V#,a)) = exp - (|a(t)| 2 - \a(0)\ 2 + a 2 (t) - a 2 (0)) +7 / ds a(s)Re (Fe~ i5s ) 



Jo 



■up 



\a{t)) 



with 



a(t) = a 



e 2* + 



^ 2i& 



Y 2i8 



idt 



As any initial state can be expanded in coherent states, with f(a) = - (a\i/;^(0)) , 



|tk(0)>= d 2 af(a)\a) 



and due to the linearity of eq. (|3.20|) , the general solution is simply 



|^(t))= j d 2 af(a)Mt,a)). 

For a single coherent state as initial condition the measured quadrature 

(^(t,a)\X\^{t,a)) 1 



(*(*)>< 



a(t)e- lip + c. c. 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



(iPt(t,a)\iPt(t,a)) 2 

turns out to be independent of the measurement noise. This is no longer true if a superpo- 
sition of coherent states is present initially. 

Using the present method the decay of the initial superposition |^(0)) = const. (\a) + 
| — a)) is studied in |2~2"]. See also pB] , [24| where this problem is studied using the nonlinear 
stochastic wave-equation. 

Finally, we determine the measure of the output process. Using eq. ( &2% with \ifj e {t)) 
given by eq. ( f4.7|) where £ is replaced by 8, one obtains a rather complicated expression. 
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Unless /(a) is an exponential linear, or quadratic plus linear, in the real- and imaginary 
parts of a, the resulting measure will be non-Gaussian. However, if (X(t))g is known as 
a functional of 6, as for the example (f4.8|) where it is even independent of 9, one may use 
eq. ( |3.33| ) to obtain 

dfi t ({e}) = dtf({0}) e*P (2V7 J\x(s)} e d6(s) - 2 1 J\x{s)) 2 e ds^ . (4.9) 

For the example ( |4.8|) the result (fOl ) is just a shifted Wiener measure. 

As a second example we turn to subharmonic generation from the vacuum, where 

H (a,a j ) = ^(a 2 -a) 2 ) (4.10) 

and 

|^(0)) = |0>. (4.11) 

The ansatz 

\^(t)) = exp (g(t) + a(t)a* + \(3{t)a ]2 ) |0) (4.12) 
solves eq. ( |3.20| ), provided the differential equations 



dj3{t) = -j/3(t)dt + (k - (k + 7e" 2 ^)/? 2 (t)) dt 

da{t) = -(l + ( K + 7 e" 2 ^)/?(t)) oc{t)dt + ^{3(t)e- lv d£(t) 

dg(t) = ~ (a 2 {t) + /3(t)) dt - ^a 2 (t)e- 2i ^dt + v ^a(t)e"^(t) (4.13) 

are satisfied with the initial conditions a(0) = (3(0) = g(0) = 0. The integration of eqs. fl4.13|) 
yields 



= K + (-I + vt + k2 + K7e " 2 ^ tanh V T + k2 + ^ e ~ 2 ^ - *o) 



a(*) = ^7 re^ (M " (K+r " 2!V)/ ^ (s) ^(r)e^^(r) 

JO 

1 /■* ft 



9(t) = -- I rfr («/3(r) + (k + 7 e- 2i > 2 (r)) + ^7 ^ a(r)e^ d£(r). (4.14) 
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The coefficient (3{t) describes squeezing. It contains the constant of integration to, which is 
fixed by the initial condition (3(0) = 0. The squeezing coefficient (3(t) is influenced by the 
measurement through 7, but unlike the coefficient a(t) which determines the amplitude and 
is completely generated from the output-noise, the squeezing coefficient f3{t) is independent 
of the noise. 

It is easy to calculate the prediction for the measured quadrature (X(t))^ as a functional 
of the output noise, here still written as £, for simplicity. We obtain 

/Y( ,^ Re[(a(t)+(3(t)a*(t))e-^] 

where a is a complex linear functional of £. In fact (X(t))^ is also a linear functional of £. 
The measure for the output noise is given by eq. ([3.271) an d may be written as 

dfH{{(}) = dfx t w m)e 9(t)+9 * (t)+la{t)l \a^ (4.16) 

The expectation value with respect to coherent states can be evaluated explicitely by the 
general form of the Baker Haussdorff formula, but we will not go into the necessary algebra 
here. In fact, like in the first example there is an easier way directly based on eq. fl3.33| ) 
and the formula ( |4.15[ ), directly yielding the Gaussian measure ( |4.9[ ) with (X(t))g given by 
( |4.15| ), after £ is replaced by 9. 

V. HETERODYNE MEASUREMENTS 

In heterodyne measurements the frequency Ui Q of the local oscillator differs from that 



of the quantum oscillator, Q = lu — u>\ 7^ 0. Therefore the phase tp in eq. ( |3.6| ) must be 
replaced by ip + fit. In particular we now define 



dG(t) = ^dB mt (t)e-^- iQt + dBt nt (t)e^ +iUt 
dE(t) = (dB- m (t)e- i<f - mt + dBljt)e itp+mt 



X(t) = - ( ae-^- im + o)e^ +iUt ) (5.1) 



18 



and eqs. ( |3.7| ), (|3.12|) hold unchanged. The entire derivation of the preceding section leading 



up to eq. ( |3.20| ) may then be taken over, and we obtain 

d\^(t)} = - iH (a, a\ t)dt - ~aW* + y/^ae' iip - int d^(t) \ip^t)) (5.2) 
which differs from eq. ( |3.20| ) only by the replacement 

d£(t) -> d£{t)e~ int . (5.3) 

The examples of section 4 can be extended to the case of heterodyning. However, it is 
not enough to make the replacement Q5.3|) in the solutions, because the Ito-rule (<i£(t)) 2 = dt 
is not analytic in d£(t). 

The wave-function ( fOj) of the harmonic oscillator is extended to the case of heterodyning 
by the replacement tp — > <p(s) = p + Qs. As the measured quadrature is independent of the 
noise it changes only because of the changed definition (|5.1|) , again leading to the replacement 
cp -> if + Vlt in eq. (f£gp. 

In the example of subharmonic generation, the generalization to the case of heterodyning 
is less straightforward. In eqs. ( |4.13| ) it leads to the same replacements (p — > (p+Qt as before. 



As a consequence the equation for the squeezing parameter (3(t) can no longer be integrated 
by a separation of variables. However, while j3(t) is changed it still remains unaffected by 
the measurement noise. Furthermore, the solutions for a(t) and g(t) in eq. (|4.14j) remain 
valid if the change <p — > p + Qt is made there. In the case where \Q\ is large compared to 
the characteristic inverse time-scales of the system it is possible to average over many cycles 
of Q 0,0]. Then 

is a complex white noise with real- and imaginary parts d£i, d^ 2 satisfying 

<%1 = d£l = \dt, di x di 2 = 0. (5.5) 



The dependence on p and f2 in eq. ( p.2[ ) then disappears. The two examples of section 4 then 
simplify again: In the oscillator wave-function (|4.4j) the term a 2 (s)e~ 2lip<yS ' > averages out and 
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disappears, and d£(s)e tip ^ = d£i(s) + id£ 2 (s) is now complex white noise. The measured 
quadrature (|4.8|) is unaffected by these simplifications and produces a signal at frequencies 
±0. 

In the example of subharmonic generation all terms with e~ 2lip ^ in eq. ( |4.13|) disappear 
by averaging out and e~ llfi ^d^(s) becomes complex white noise. The solution is then given 
by (|4.14|) without the e~ 2iv -terms and with complex noise. The phase (p(t) — tp + Qt then 



only remains in (X(t))^ in eq. ( 4.15 ) and again produces a signal at ±f2. 



VI. PHOTON COUNTING 

In this section we consider the case of photon counting by an ideal photodetector placed 
behind the partially transmitting mirror of a cavity. The concept of input and output 
processes used in the preceding section to describe homodyne and heterodyne measurements 
can be generalized to this case. Let us introduce the input photon number process A in (i) 
defined by (see e. g. [24-26]) 

Ajn(t) = f bl(s)b in (s)ds (6.1) 
Jo 

where formally 

dB in (s) 

Oin(s) = — j • (6.2) 

as 

The increment <iA in (t) can be expressed in terms of the white-noise creation and annihilation 
operators as 

dAln(( ) = ^y-W, (6 .3) 
It satisfies the commutation relations 

[dA in {t),dA in {t')\ = 0, [dA in {t),X(s)\ = for t > s (6.4) 
for any system variable X. We shall denote the eigenvalues of c?A in (t) by d\(t). 
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The output number process can be calculated using the Hamiltonian H mt in the interac- 
tion representation of eq. ( |3.3| ) as a unitary transformation 

Aout(t) = uUt)Mn(t)U int (t) (6.5) 

of the input number process, with U- mt defined in eq. Q3.4p. The increment cfA out of the 
output process satisfies the quantum stochastic differential equation [24-26] 

dA ont (t) = dA- m (t) + ^{a{t)dB\jt) + dB- m {t)a\t)) + ja\t)a(t)dt (6.6) 

where all system operators are in the Heisenberg picture. Note that the output photon 
number process satisfies the same equal-time commutation relations and multiplication rule 
as the input process fl25| , |27jl . Eq. (|6T6|) is the analogue of eq. ( |3.7p . The difference lies in the 
fact that the output number process does not commute with the input number process for 
equal times as can be seen from eq. flBTB]) and the commutation relation [dB in (t), dBj n (t)] = 
dt; as a consequence the input and output number processes at equal time can not be 
diagonalized simultaneously, i. e. it is not possible to express the eigenvalues of the output 
process in terms of the d\(t). 

Let us consider now the time evolution of the state of the total system. The generator 
L(t)dt is unchanged and still given by eq. (|3.16|) . Up to now we have not specified the 
initial state of the bath which is necessary to evaluate the term dB\ n (t)dB- m (t) . Assuming 
the bath to be in a vacuum state as in section 3 we immediately see from eq. (|6.3| ) that 
the input number process vanishes, which is clear from a physical point of view: due to its 
absorbing nature a photodetector cannot 'see' the vacuum fluctuations. The mean of the 
output process is given by (dA ont (t)) = j(ala)(t)dt as can be seen by taking the expectation 
value of eq. ( |6.6| ) in the vacuum state of the bath. However, in this case it is not clear how 
to apply the input-output formalism to derive the equation governing the time evolution 
of the state of the system. A calculation of this time evolution can still be performed in a 
manner similar to the preceding section if we take the bath to be in a coherent state |{/?}) 
with 
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dB in (t)\{(3}) = ^e(t)dt\{(3}) (6.7) 

where yffe(t) is the arbitrary complex amplitude of the coherent state. At the end of the 
calculation we may take the limit e — > 0. Physically, this assumption means that not only the 
vacuum fluctuations but also a small classical field enter the cavity of the quantum oscillator. 
According to eq. (|6.7|) the product dB\ n {t)dB m (t) is of second order in dt and therefore can 
be set equal to zero in the following (this substitution can be justified rigorously in terms of 
quantum stochastic integration; see e. g. p8|). We then obtain the following multiplication 
table 

dB in (t)dB in (t) = = dB\ D (t)dB\ Q (t) 

dBl(t)dB in (t) = 
dB in {t)dBljt) = dt 

dB- m (t)dt = = dB\ n {t)dt 

dA in {t)dB in {t) = 
dA- m (t)dA in (t) = dA- m (t) 

dA in {t)dt = 0. (6.8) 

For the eigenvalues of dA in (t) = d\(t) the condition 

{d\{t)) 2 = dX(t) (6.9) 

implies 

d\(t) = < . (6.10) 

I 1 

The dX(t) at different times are statistically independent due to the statistical independence 
of the underlying quantum Poisson process at different times. 

Proceeding as in section 3 we compute the wave function of the total system at time t as 

t 



\Aot(t)) 



H (1 + L(U)dti 



|V(0)> ® |{/9}> (6.11) 
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where we assume that initially the state of the total system factorizes. Now the operators 
dB- m (ti) can be commuted to the right where they are replaced by ^e(ti)dU when acting 
on | {/?}). Then the Schrodinger equation for the total wave function is given by 

\d + iH Q {a,a\t)db + |aWt - ^(adB} n (t) - ^ye(t)aUt)^j \if>(t)) ® \{/3}) = 0. (6.12) 

Now we make the transformations 

1 adA- m (t)m)) ® \{(3}) (6.13) 



e(t) 

where we have used eq. ( |6.3| ) in the last step. In the representation where the input number 
process dA. in (t) can be replaced by the eigenvalue d\(t) the Schrodinger equation flS.12| ) is 
therefore given by 

t)dt + -a ] adt + yefyaUt - —ad\(t) ) \ip\(t)){\\{(3}) = (6.14) 
\ 2 e(t) J 

so that the bath state can be divided out. In this way we obtain a linear stochastic 
Schrodinger equation for the state of the system alone 

d\ip x {t)) = \-iH (a,a\t)dt - ^aWt - ye(t)a^dt + —ad\(t)\ \if>\(t)) (6.15) 
I 2 e(t) J 

with the measure 

^ P ({A}) = |(A|{/3})| 2 . (6.16) 

Here and in the following |A) denotes the number states of the bath satisfying 

dA in (s)|A) = dA(s)|A) (6.17) 

for all s. Discretizing time t — *■ ij and using the fact that d\(ti) is or 1 we obtain the 
formal expression 

N-l 

dtf ({A}) = J] [C 1 - 7|e(ti)| 2 ^)(l - dX(U)) + ^W^dtdKU)] . (6.18) 

i 
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Thus d\(t) is an independent Poisson process with mean value 

(dX(t)) = ({{3}\dA in (t)\{[3}) = i\e(t)\ 2 dt. (6.19) 

Using eq. ( |6.15| ) and the fact that d\(t)dt = we can write the wave function of the system 
at time t + dt in the form 



\Mt + dt)) = \Mt)) + d\Mt)) 



= (1 - dX(t)) (1 - iH (a } a\ t)dt - -^a)adt - ^e(t)a ] dtj \^\{t)) 

^A(f)^ + lW(f)>. (6.20) 

This expression clearly shows how the wave function evolves during the time interval dt 
depending on the result of the photon counting: if no photon is detected (with probability 
1 — 7|e(t)| 2 <i£ according to eq. ( |6.18|) ), i. e. dX{t) = 0, this time evolution, according to 
eq. ( |6.20D , is governed by the non-hermitian Hamiltonian H = Hq — i^a^a — i"je(t)w\ if a 
photon is detected (with probability 7|e(t)| 2 <it), i. e. dX(t) = 1, the wave function jumps 
according to the second term in eq. ( |6.20p . This time evolution of the wave function due to 
the linear stochastic Schrodinger equation fl6.15| ) is analogous to the nonlinear prescription 



introduced phenomenologically by Carmichael |15| to generate numerically what he has 
called 'quantum trajectories' for the wave function of the system conditioned on the measured 
photon number (see also |J). Applying eq. ( |6.20| ) repeatedly in time the formal solution of 
the linear stochastic Schrodinger equation (|6.15| ) given by eq. ( 6.1 1| ) can be written in the 
more elegant form, with discretized time t i: t = 0, tjy = t 
ryv-i 



\Mt)) 



i=0 

+dX(U) (-?— + 1 



JJ "I (1 — dX(ti)) yl — iHo(a, a\ t)dt — —a^adti — ^t{ti)a) dt,, 



|V(o)> 



= K x (t, 0)1-0(0)) (6.21) 

which will be useful for the calculation of correlation functions. It is important to note that as 
a consequence of factoring out the state of the bath in eq. (|6.14|) the time-evolution operator 
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K\(t, 0) is non-unitary. It is now straightforward to obtain the analogue of eq. (|3.23|) for 
arbitrary multi-time correlation functions with dfi w replaced by dfi p . 

Let us show now how the complete information on the output process is contained in 
the solution of eq. ( |6.15| ). To this end we define a jump process n(t), with (d,K,(t)) 2 = ck(i), 
which is distributed according to the measure 



dtl t ({K})=d^({K}){Mt)\Mt)) 



(6.22) 



where d\i P ({k}) is the Poisson measure (|6.18|) of the input process. The measure ( |6.22|) is 
normalized, i. e. 



dfH{{ K }) = I d{K}\({(3}\K)\ 2 (Mt)\Mt)) = J d{n} (Mt)\ ® ({/3}k><«|{/3}> ® \Mt)) 

d{K}(^(*)l«X«lV'tatW) = (Aot(t)\Aot(t)) = 1 (6.23) 



because the total wave-function of system and environment is normalized. Going through 
the same argument as in section 3, eq. fl3.27|) - fl335"D , we obtain for the conditional measure 
of a single time step 



p(dK(t)) = p F (dn(t)) 1 + 



(Mt)\Mt)) 



(6.24) 



where 



p p {dK{t)) = (1 - 7|e| 2 dt)(l - dn(t)) + ^\e\ 2 dtdn{t) . 



(6.25) 



The increment of the norm of \ip K (t)) can be calculated from eq. ( 6.20|) and we obtain 



p(dK(t) = 0) = (1 - -i\e(t)\ 2 dt) (l - 7(e(t)a t + e*[t)a) K dt - ~{{a)a) K dt 
= 1 - j\e(t)\ 2 dt - 7(e(t)a f + e*(t)a) K dt - -f(a j a) K dt 



p(d«(t) = 1) = iW)\ l dt 1 + 



+ 



+ 



v <t) e*(t) \e(t)\ 2 
j\e(t)\ 2 dt + 7(e(t)a f + e*(t)a) K dt + -y(a^a) K dt 



- 7(e(t)a t + e*(t)a) K dt - ^f(a^a} K dtj 

(6.26) 



where we defined for an arbitrary operator Q 
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m)) _ (Mmm (6 27) 

The expectation value of dn{t) according to this normalized probability distribution is simply 
given by 

(dn(t)) = p{dn{t) = 1). (6.28) 



Comparing eq. (|6.26 ) with the expectation value of eq. ( |6.6| ) in the coherent state of the 



bath we see that the expectation value of the process dn(t) introduced by eq. Q6.22 ) is equal 
to that of the output process dA ont (t). 

We can go even further and show the equality of all correlation functions of dn(t) and 
the corresponding correlation functions of dA out (t), which means that dn(t) is stochasti- 
cally equivalent to the output process. The calculation of correlation functions of dn(t) is 
greatly simplified by the use of the formal solution ( |6.21| ) of the linear stochastic Schrodinger 



equation ( |6.15| ) (see also section 3). Let us consider first the two-time correlation function 
/ dfi t+T ({K})(dK(t + r)dK{t)) = (d«(t + r)ck(*)) (6.29) 



with r > 0. According to eq. ( |6.22| ) we can reexpress eq. ( |6.29 ) in terms of the two-time 



correlation function of the input process as 

(dK{t + r)dre(t)) = J dfi( +T {{X}) d\{t + T)d\(t)(ip x (t + r + rfr)|^ A (t + r + dr)) (6.30) 



where \ip\(t + t + dr)) is given by eq. (|6.21 ). Now we assume that the wave function at time 



t is given and not conditioned on the input process, so that \if)(t)) is simply the initial state 
of the system. The components of \ip\{t + r + dr)) with d\(t) = or d\(t + r) =0 cannot 
contribute in the integral ( |6.30| ). Therefore we can rewrite eq. ( |6.30| ) as 



(dK(t + T)dK{t)) = J dtf +T {{X}\dX{t) = dX(t + r) = 1) (ip x (t + r + dr)\ip x (t + r + dr)) x 
xp p {dX{t) = l)p p {dX{t + r) = 1) (6.31) 

with the Poisson measure conditioned on dX(t) = dX(t + r) = 1 and 
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\Mt + r + dr)) = + 1 j Kx(t + r,t + dt) + lj (6-32) 

The non-unitary evolution operator K\(t + r, t + dt) is defined in eq. (|6.21| ) and we have 
used again the fact that d\(t) = d\(t + r) = 1. Since the probability p p (d\(t) = 1) is 
independent of the state \ip\(t)) of the system we obtain 

(dn(t + T)dK(t)) = 1 2 \e{t)\ 2 \e{t + r)\ 2 dt 2 J dtf +T ({V)(m\ (j^ + l)Kl(t + r,t + dt) X 



+ 1 ) (^r^ + ^Kx(t + T,t + dt)(^- + i)m)) 



\e*(t + r) J \e{t + r) J M ' \e(t) 

= 1 2 dt 2 J drf +T ({\}) (ip(t)\(a^ + e\t))K\it + r,t + dt){<J + e*(t + r)) x 

x(a + e{t + r))K x (t +r,t + dt)\ip(t)). (6.33) 

Due to the non-unitarity of K\(t + r, t + <it) it is not possible, in general, to insert K\K\ 
between factors at will. However, because the total time evolution of system and reservoir 
is unitary, we still have unitarity under the integral 

/ drf +T ({\}) K{{t + r,t + dt)K x (t + r,t + dt) = l (6.34) 

as may be checked by using eq. ( |6.21|) . Therefore under the integral over dfi p ({\}) the 



stochastic time-evolution operator K\(t + r, t + dt) may be treated as unitary, and in this 
sense K x {t) can be used to define a stochastic Heisenberg picture, 

Q A (t) = K\(t, + dt)QK x (t, + dt). (6.35) 

Using this we may now rewrite eq. ( |6.33j ) as 

{dn{t + T)dn{t)) = ^ 2 dt 2 f dfif +T ({X}) (ip\a\(t)a{(t + r)a x (t + r)a A (t)|^) (6.36) 



where we used a(t) = a(t) + e(t). The right-hand side is the standard form for the degree 
of second order coherence g^(t,t + r) to be found in textbooks on quantum optics (see 



e. g. @). 

The above procedure can be generalized in a straightforward manner for the calculation 
of correlation functions of higher order. The general n-time correlation function for unequal 
times ti tj for all i, j is thereby obtained as 
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(dK(t n )...d K (t l ))= 1 n dt n I dy^ n {{\}) (i/)\a\(ti) ■ ■ ■ a\(t n )dx(t n ) ■ ■ ■ d\(ti)\ip) (6.37) 



where we have assumed that t n > t n _i > . . .t\. Eq. (|6.37|) is equal to the degree of n-th 
order coherence. 

To complete the comparison of the process dn(t) with the output process dA out (t) we 
compute the general normal ordered ra-time correlation function of the output process 

(: dA ont {t n ) . ..dA ont {h) :) = (^i)l <g> ({/3}| : dA out {t n ) . ..dA ont (h) : \{p}) <g> (6.38) 

where : : denotes normal and time-ordering. Note that the expectation value is taken in the 
state of the total system because the output process dA ont (t) acts in the Hilbert space of 
the total system. Expressing dA out (t) in terms of the input field by eq. ( |6.6D and using the 
statistical independence of dA- m (t), dB in (t) and dB\ n (t) at different times we finally obtain 

(: dA ont (t n ) . . . dA out (h) :) = 7 n rft n (a t (t 1 ) . . . a\t n )a(t n ) . . . a(*i)) (6.39) 



where the expectation value on the right side of eq. ( |6.39 ) is again taken in the state of the 



total system. Now comparing eq. ( |6.39| ) with eq. (|6.37|) we see that both expressions are 
equal and dn(t) is indeed statistically equivalent to the output process. 

Let us now make contact to the nonlinear stochastic Schrodinger equation used by others. 
This is done by replacing the input process d\{t) in eq. ( |6.15| ) by the output process dnit), 
i. e. 

d\{t) dn{t) _ dn(t) 



yj{d\{t)) ^{l{W)\ 2 + (e(t)at + e*(t)a) K + (rta) K ))dt y/P(t)dt 
where P(t)dt = (dn(t)). This replacement is phenomenological in distinction to the micro- 
scopic derivation of eq. ( |6.15| ) we have given above. It corresponds to the transition from 
the linear equation ( |2.2| ) to the nonlinear equation ( |1.3| ) by taking as equal to the ouput 
probability \(ijj^\ijj)\ 2 . In this way we obtain a nonlinear stochastic Schrodinger equation for 
the non- normalized wave function \ip K (t)) 

d\ip K (t)} = \ -iH (a, a\ t)dt - -ctadt - r ye{t)a s dt + Jj ,° dn(t) \ \i/j K (t)). (6.41) 
I 2 \fW) J 
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The nonlinearity arises because P(t) depends on the wave-function due to the appearance of 
expectation values in its definition ( |6.40| ). The superposition principle is thereby destroyed. 
Performing the limit e(t) — > 0, i. e. setting the coherent driving field equal to zero this 
stochastic Schrodinger equation takes on the form 

d\il> K (t)) = \ -iH (a, a\ t)dt - ^a ] adt + = dn(t) \ \if) K {t)) (6.42) 



and for the normalized \<j> K (t)) = IM*)) / (MWMW 



d\(f> K (t)) = { -iH (a, a\ t)dt - -(a f a - (a f a))dt + I = - 1 I ck(f) ) \4> K (t)) (6.43) 



2 



(ata) K 



which is identical to the nonlinear stochastic Schrodinger equation introduced by Carmichael 
| 15| (see also [p]]). Note that the classical output process cannot be reexpressed in terms of 



the classical input process as was already stated in the beginning of this section. It should 
also be noted that the limit e — > obviously cannot be taken in the linear equation ( |6.15D . 



However, even in eq. (|6.41|) this limit cannot be taken in a strict form, because the resulting 
eqs. ( |6.42| ), ( |6.43| ) are ill defined for the vacuum state, where (a^a) K = 0. 



Finally, for illustration we apply the linear stochastic Schrodinger equation ( |6.15| ) to the 
periodically driven harmonic oscillator already studied in section 4, i. e. we solve the linear 
stochastic Schrodinger equation 



j_2 (F*e i5t a - Fe~ iM a ] ) dt - ^adt - 1 e(t)aUt + -^-d\(t)\ |^ A (t)) (6.44) 
[2 2 e(t) J 



where initially \ip(0)) = \ao). Inserting the ansatz \ip\(t)) = g(t)\a(t)) in eq. (|6.44|) gives the 
following two differential equations 



dg(t) 



-~ {\a{t)\ 2 + ilm (F^a(t)e i5 ') + 2Re(a*(t)e(t))} dt + 

2 €{t) 



g(t) (6.45) 



with initial conditions g(0) = 1, a(0) = ao- The integration of the differential equation for 
a(t) yields 
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,,(/) ^ ! " ifw ) + I3W e ^ -ie~^ Jldse{s)ei s (6.46) 

7 / 7 ' ' 

while the solution for can be written in the manner of eq. ( |6.21| ) as 

AT-l 

n 

i=0 



9(t) = II (l~dA(ti)) f-|{|a(^)| 2 + zIm(F*a(t i )e i ^)+2Re(a*(^)e(t i ))}^ 



with t = 0, £at = £. 



(6.47) 



VII. CONCLUSION 

Traditionally the time evolution of quantum systems comes in two forms, one is linear, 
unitary, and deterministic, and describes systems which are completely isolated, even from 
any external measuring apparatus; the other is nonlinear, non-unitary and random and 
describes the results of measurements. In the present paper we have discussed how the non- 
linearity in the description of measurements can be avoided, first for the usual schematic and 
idealized general description of quantum measurements, and then more concretely for the 
well-known quantum optical measurement schemes of homodyne- and heterodyne measure- 
ments and photon-counting. In each of these concrete examples we have given microscopic 
derivations of the linear wave-equations, which govern the measured systems. In the limit 
of vanishing coupling with the measurement devices (7 — » 0) they reduce to the Schrodinger 
equation for closed systems. The existence of these equations should dispel the old believe 
that the two forms of quantum dynamics mentioned above are incompatible. Rather they 
appear here as limiting forms of the microscopically founded description we have given. This 
is only possible because a linear description of non-isolated measured quantum systems ex- 
ists, which preserves the validity of the superposition principle also for this class of systems. 
Besides some practical advantages in certain special cases this is the fundamental reason 
why we have put all emphasis on the linear form of the wave-equations and mentioned their 
equivalent nonlinear counterparts only in passing. In this linear description the effects of 
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the environment or the measurement apparatus resides only in the non-unitary dissipative 
terms and the non-unitary stochastic terms of the wave-equation. These alone are sufficient 
to produce the familiar disappearance of interference terms as the coupling of the system to 
the external world is increased, a nonlinearity of any form is not required to produce this 
effect. 

For the three measurement schemes which we have discussed in detail we have shown 
that the solutions of the linear wave-equations yield 

• expectation values and correlation functions of all system variables from formulas 
like eq. ( |3.23|) , where the unnormalized stochastic expectation value taken with the 



stochastic wave-function is averaged over the Wiener noise, in the case of homodyne 
or heterodyne measurements, or over the Poisson noise in the case of photon counting; 

stochastic realizations according to their correct measure of the measured output field 
0(i) given by eqs. ( |3.6| ) or (|5.1| ), or A out (t) of eqs. ( |6.5| ), flBUD in the case of photon- 
counting, via the normalized stochastic expectation values (X(t))o or (a' (t)a(t)) K in 
eqs. ( |3.33| ) or ( |6.26| ), respectively; this result is particularly important because it allows 



to evaluate all correlation-functions, power spectra, etc. of the measured quantities 
from a single, sufficiently long time series via eqs. (|3.37| ) or (|6.3?1 ), respectively. 



• the complete measure of the output field of the measurement via eq. (|3.27|) , or 
eq. (|6.22| ), respectively. 

The practical advantage of the linear description we have illustrated by giving exact so- 
lutions for a number of examples. Of course many more examples for exact solutions can 
and undoubted will now be given: all examples of nonlinear quantum optics which allow 
for a linearization (e. g. linear parametric amplifiers, Raman amplifiers, nonlinear quantum 
oscillators like lasers if operating far above threshold) fall in this class. However, the linear 
stochastic wave-equations present crucial advantages even if intrinsicly nonlinear physical 
processes are considered whose quantum effects under measurements or interactions with 
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reservoirs are of interest — e. g. the important class of nonlinear mesoscopic systems exhibit- 
ing at the same time quantum interference effects like dynamical localization or Ericson-type 
fluctuations and traces of classical chaos. Obviously, if stochastic wave-equations are chosen 
to describe such systems it is only their linear version which will permit to give a transparent 
description of such quantum interference effects. For instance path-integral solutions of the 
linear wave-equations exist and semi-classical approximations can be based on them. 

So far the theory has only been developed within the Markovian framework. In fact, 
the reliance on the Markovian limit is heavy and one may suspect that this limit is really 
essential in order to achieve the description of the subsystem by a stochastic wave-function. 
Although this may be possible, it has not been shown, and non-Markovian generalizations 
of the formalism we have presented remain a challenge. 
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